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Abstract 

Traversal of a symmetry-breaking phase transition at a finite rate can lead to causally- 
separated regions with incompatible symmetries and the formation of defects at their bound- 
aries. The defect formation follows universal scaling laws prescribed by the Kibble-Zurek 
mechanism (KZM) important to the study of phase transitions in fields as diverse as quan- 
tum and statistical mechanics, condensed matter physics and cosmology. Here, we observe 
the KZM in a crystal of cold trapped ions, which is conducive to the precise control of 
structural phases and the detection of defects. The experiment confirms a scaling law with 
an exponent of 2.68 ± 0.06, as predicted from the KZM in the finite inhomogeneous case. 
Such precision makes it feasible to use ion crystals for quantitative tests of classical and 
quantum statistical mechanics. 
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T. W. B. Kibble 01 [2 proposed to model the emergence of domain structures in the first 
moments of the cosmic evolution by a symmetry-breaking phase transition. Later, W. H. Zurek 
suggested laboratory tests of this mechanism through the observation of a quench-induced phase 
transition in condensed matter (3). Experiments using liquid crystals flU and liquid 4 He (H 
|6l 13 provided the first observations of the KZM. Recent theoretical studies jH HI [TOj have 
triggered further interest in studying KZM scaling laws in tailored, well-controlled systems 
with high precision. Here, we present an experimental observation of the KZM, using a laser- 
cooled Coulomb crystal of 16 calcium ions in a linear Paul trap. The simplicity of this finite 
system permits the observation of the KZM under very clear experimental conditions without 
disturbances from initial defects or impurities. Moreover, the experiment features a highly 
controlled environment with the possibility to tune the control parameters over a wide range, 
ideal for the study of non-equilibrium statistical mechanics. 
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Figure 1: Fluorescence images of a 16 ion crystal during the measurement of the KZM. a) Linear 
ion crystal before ramping the axial potential; b), c) zig-zag/zag-zig configuration after the ramp; d, e) 
appearance of single defects, which connect incompatible orientations of the crystal; f) double defects 
within the crystalline structure. The red line clarifies the configuration of the crystals. The exposure time 
for each of the images was 30 ms. The width of one pixel corresponds to 2.2 fim. 

The KZM applies to non-equilibrium systems traversing a second-order phase transition. 
Prior to the phase transition, perturbations relax into the lowest energy equilibrium state under 
the dissipative influence of a cooling mechanism. This defines a characteristic relaxation time. 
Near a structural phase transition, this relaxation time increases as it diverges at the critical point 
(CP). Therefore, there is no finite rate at which the CP could be traversed adiabatically and the 
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structure of the system is effectively frozen even before it reaches the CP. Furthermore, different 
regions of the system are causally disconnected due to a finite propagation of information, 
allowing for the presence of different symmetry-broken ground states. If the choice of symmetry 
of neighboring sections are incompatible, defects form where the phase boundaries meet and 
the system is thus prevented from reaching a global ground state. 

For an inhomogeneous system, different regions reach the CP at different moments and 
lead to time dependent phase boundaries. This gives rise to an adiabatic transition if the speed 
of information transfer v s exceeds the speed of the propagation of the phase boundary v p , thus 
allowing for the causal connection of different domains and preventing the formation of defects. 
Hence, defects may only be created in a region of the system where v p > v s . 

Laser-cooled Coulomb crystals in ion traps represent finite-sized systems with a tunable 
spatial inhomogeneity. We investigate the occurrence of the KZM in the structural phase transi- 
tion from the linear to zig-zag configuration l03l[T4ll . The regular structure of localized ions in 
a linear Paul trap is controlled by the confining electrostatic potentials along the length of the 
crystal (axial) and the dynamic radio-frequency potential in the perpendicular (radial) directions 
which gives rise to an effective harmonic potential. The combination of the Coulomb repulsion 




Figure 2: An X-shaped micro-fabricated segmented Paul trap based on four laser-cut alumina chips 
each with 1 1 electrodes. The radial confinement is generated by applying a radio-frequency voltage of 
Urf ~ 450Vp P at a drive frequency of Q/(2ir) = 22 MHz, resulting in a relevant radial trap frequency of 
w ra d/ (27r)= 1.4 MHz. The axial potential is generated by a superposition of static DC potentials applied 
to the segmented electrodes and variable voltages applied to the conical end-cap electrodes allowing for 
axial frequencies ranging from w ax /(2vr) = 160kHz to 350kHz (see Fig. ^. A field -programmable- 
gate-array (FPGA) enables precise timing control for triggering the camera system and an arbitrary 
waveform generator, which has 16 bit amplitude resolution and 100 MS/s for fast and glitch- free voltage 
ramps with variable time constants lfTTl[T2l . Optical detection is achieved using an electron multiplying 
charge-coupled device camera with a 10-ms-exposure time, oriented at 45° to the planar structure of 
the crystal. The blue arrow indicates the direction of the Doppler cooling laser, with projections in all 
principal axes of the trap. 
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between ions and the dissipative forces of laser-cooling yield a regular crystalline structure with 
minimum potential energy. For low axial trapping frequencies, this structure is a linear chain 
(Fig. [T]a). When the strength of the axial confinement is increased beyond a critical value, the 
structure is squeezed into the radial dimensions and undergoes a structural phase transition from 
a linear crystal to a zig-zag crystal, which can be described by a second-order phase transition 
jH. The latter phase can take the form of two symmetry-broken ground states, referred to as 
"zig-zag" (Fig.[Tj)) and "zag-zig" (Fig.[l]:) configurations. Here the phase transition initiates at 
the center where the charge density is highest and the phase boundary propagates outwards to 
both ends with finite propagation speed v p . If the curvature of the trapping potential is rapidly 
increased, v p can exceed the speed of sound v s , and defects may be created lfT5l [T6l . Defects 
manifest as kinks in the zig-zag structure of the crystal (Fig. [TJi-f). The positions of the in- 
dividual ions are determined by fluorescence imaging, revealing the final configuration of the 
crystal. This allows for the determination of the defect formation rate d as a function of the 
rate at which the CP is traversed. By solving the time-dependent Ginzburg-Landau equation 
in the underdamped regime ifTTll . a universal power law for the density of defects with respect 
to the rate of change of the control parameter - the derivative of the trap frequency at the CP 
7 := du>/dt\ c - can be derived. If the lengths of causally connected regions are small compared 
to the system size, a scaling of d oc 7 4//3 is predicted B9KT01. However, for small systems the 
size of these causally connected regions are comparable to the crystal length and the scaling 
becomes d oc 7 8 ^ 3 , where the exponent has been doubled due to finite size effects ifTHl . 

In order to determine the scaling law of the defect formation rate, we load sixteen 40 Ca 
ions into a linear segmented Paul trap (see Fig. [2]). Doppler cooling close to saturation is ap- 
plied during the whole experimental cycle. We ramp the axial trap frequency from 167kHz to 
344 kHz (see Fig. [3]) across the CP at 202 kHz. The voltage ramp is driven at different rates 
with time constants ranging from 0.5 fxs to 4.0 /is. An image of the crystal configuration is cap- 
tured 100 fj,s after the ramp with a 10 ms exposure time (see Fig. [T]). The number of defects can 
be seen directly by manual inspection of the images. However, to guarantee an efficient and 
reliable categorization of the images into different classes of possible crystal structures, their 
normalized two-dimensional Fourier transforms are calculated. The resulting spectra are com- 
pared to those of reference images via the sum of squared residuals. The reference images were 
generated by averaging manually selected samples for a total of 14 relevant configurations. A 
recognition threshold is applied to the sum of squared residuals to reject low-quality images 
(less than 5%). The measured defect formation rate as a function of the rate of change of the 
axial trapping frequency 7 is presented in (Fig. [4]) for 60000 experiments. The predicted scaling 
law for the inhomogeneous finite-sized system is confirmed by the observation of a slope of 
(3 = 2.68 ± 0.06, matching the theoretical prediction of 8/3 ~ 2.67 over the full dynamic range 
of defect formation rates d. 

The finding is supported by a molecular dynamics simulation to analyze the ion trajectories 
during a ramp. A realistic model of the trap is employed, and the equations of motion are solved 
with a partitioned Runge-Kutta integrator |fi~9l . Laser-cooling is implemented using a Monte- 
Carlo simulation and constant dissipation. The micromotion for each particle is fully modeled 
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Figure 3: Top: Measured voltage applied to the end-caps by the arbitrary waveform generator (black 
line). A shape of the functional form V(t) oc [1 + exp(— (t — io)/ T ))] 1 is chosen (red line) to avoid 
excitation of axial vibrations. The timescale parameter r determines the rate of change of the con- 
trol parameter 7 at the CP. The dashed line shows the separation between the two structural phases at 
w a x/(2vr) = 201.7 kHz and a radial trap frequency of w ra d/(27r) = 1394.1 kHz. Inset: Dependency of 
the trap frequency / ax on the applied end-cap voltage. A square-root function fits the measured data. 
From this measurement, the functional dependency of the trap frequency on time and thus the rate of 
change of the axial frequency at the CP is deduced. Bottom: Axial positions of the ions during the ramp 
as extracted from simulation results. Diamonds indicate the onset of the local phase transition for each 
ion, which are reached at different times due to the inhomogeneous charge density. The dashed line 
indicates the time when the middle ions reach the CP. After the compression, the motional energy of the 
crystal is slightly increased, however, not sufficiently to cause the ions to swap positions. 

by including the oscillatory electrode potentials [|20l . By monitoring the trajectory of each 
individual ion, we verified that no swapping of ions, and therefore no melting of the crystal, 
occurs during the ramping procedure and ensured that the excitation of axial oscillations is 
minimized (see Fig. [3]). 

In conclusion, we have observed the KZM in a model system with ideal preparation, control 
and readout capabilities. The observed scaling of the defect formation rate is in excellent agree- 
ment with the theoretical prediction. In future experiments, we might utilize the trap control 
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Figure 4: Double-logarithmic plot of the measured rate of defects d over the rate of change of the 
axial trap frequency 7 = (dui ax /dt) | c at the critical point. The power law fit of d oc 7^ (red line) 
results in an exponent of /3 = 2.68 ± 0.06, which is in agreement with the theoretical prediction of 
/3 = 8/3 « 2.67. Only data points (l-cr error bars) within the solid red line are used for the fit. The 
uncertainty in (dtu a x/dt) \ c is deduced from the scatter of repeated recordings of voltage ramps. The 
constant offset visible at lower ramping rates stems from background gas collisions at a base pressure of 
1 x 10~ 9 mbar. This is supported by the observation that the background rate is directly proportional 
to the measurement time. The saturation is due to the maximum number of defects being limited in the 
finite system. 

voltages to modify the local charge density and explore the role of spatial inhomogeneities for 
defect formation. Instead of external trap potential ramps, spin dependent forces could be used 
to initiate structural phase transitions and quantum quenches fl2Q. Ultimately, it might be possi- 
ble to cool the system deeper into the quantum regime to explore quantum statistical mechanics 
where phase transitions are driven by quantum rather than thermal fluctuations E2ll23ll . 
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